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Using state-of-the-art many-body calculations based on the L GW plus cumulant’ approach, we 
show that electron-plasmon interactions lead to the emergence of plasmonic polaron bands in the 
band structures of common semiconductors. Using silicon and group IV transition metal dichalco- 
genide monolayers (AX 2 with A=Mo,W and X=S,Se) as prototypical examples, we demonstrate 
that these new bands are a general feature of systems characterized by well defined plasmon reso¬ 
nances. We find that the energy vs. momentum dispersion relations of these plasmonic structures 
closely follow the standard valence bands, although they appear broadened and blueshifted by the 
plasmon energy. Based on our results we identify general criteria for observing plasmonic polaron 
bands in the angle-resolved photoelectron spectra of solids. 


Electron-boson interactions are pervasive in many- 
body physics and the resulting quasiparticles are clear ex¬ 
amples of emergent behavior in quantum matter. While 
in the case of phonons and magnons the fingerprints of 
their interactions with electrons in angle-resolved pho¬ 
toemission spectra (ARPES) are largely understood [T], 
much less is known about electron-plasmon interactions. 
A detailed and quantitative description of these inter¬ 
actions is key to refine our understanding of electronic 
excitations in condensed matter, and could provide new 
pathways towards plasmon-assisted band-gap tuning [2], 
or the manipulation of plasmon polaritons, with poten¬ 
tial implications for photonics and plasmonics [3j. 

In ARPES, the acceleration of a photo-electron upon 
photon absorption may trigger shake-up excitations in 
the sample, leading to the emission of phonons, electron- 
hole pairs, and plasmons, the latter being collective 
charge-density fluctuations [4j. Intuitively, if a photon 
excites both a hole and a plasmon, the ARPES signal 
should exhibit spectral weight at energies corresponding 
to the sum of the binding energy of the electron and the 
excitation energy of the plasmon, as obtained for example 
from electron energy loss spectroscopy (EELS) [5]. This 
phenomenology is analogous to the emergence of ‘peak- 
dip-hump’ structures in ARPES as a result of electron- 
phonon interactions OQUHIEI; the difference between the 
resulting spectral features arises from the characteris¬ 
tic energy of the boson (~10 meV for phonons, ~10 eV 
for plasmons). In fact the very first model of electron- 
plasmon interactions 0HH is formally equivalent to 
the electron-phonon Hamiltonian developed for the po¬ 
laron problem M- In this model, the electron-plasmon 
interaction results in ‘plasmonic polarons ’ m , in com¬ 
plete analogy with the polarons of the ordinary theory of 
electron-phonon interactions [16] . 

Identifying plasmonic polarons in ARPES spectra is 
notoriously difficult. While plasmonic satellites have 
been successfully identified in the integrated photoemis¬ 
sion spectra of Na m and Si [18], the identification of en¬ 
ergy vs. momentum dispersions of plasmonic polarons in 
angle-resolved spectra has proven considerably more chal¬ 


lenging [IS]. So far such dispersions have been observed 
only in the case of graphene, and only in a narrow region 
of the Brillouin zone around the Dirac point [20] EE]- Key 
factors hindering the observation of the dispersion rela¬ 
tions of plasmonic polarons are (i) the energy scale of the 
plasmon energy, which requires using energetic photons 
at the expense of momentum resolution; (ii) the increased 
phase-space for electron-phonon scattering and electron- 
hole pair generation, which adds to the spectral broad¬ 
ening; and (iii) the possible mix up of weak plasmonic 
satellites and strong quasi-particle peaks. 

In this work we perform state-of-the-art first-principles 
calculations to show that electron-plasmon interac¬ 
tions lead to the formation of plasmonic-polaron band- 
structure replica. These new structures appear as broad¬ 
ened copies of the valence bands shifted by the plasmon 
energy. Using a combination of many-body perturba¬ 
tion theory in the GW approximation [221424] and the 
cumulant expansion approach d US] HU HH1 [25H29] we 
demonstrate the presence of plasmonic-polaron bands in 
silicon. We further show that two-dimensional group 
IV transition-metal dichalcogenides (TMDs) AX 2 with 
A=Mo,W and X=S,Se [30] [31 provide an ideal play¬ 
ground for the experimental observation of these novel 
spectroscopic signatures of the electron-plasmon cou¬ 
pling. 

Within the sudden approximation, the photocurrent 
measured in ARPES experiments is proportional to the 
electron spectral function A(k, oj) pH [323 , where k is the 
crystal momentum of the electron and cu its binding en¬ 
ergy (here and in the following atomic units are under¬ 
stood) . The spectral function can be calculated by using 
the cumulant expansion d na na ns]: the electron 
Green’s function is expanded in terms of the screened 
Coulomb interaction W, and a subset of diagrams is eval¬ 
uated to all orders of perturbation d ■ This strategy leads 
to a more accurate treatment of dynamical correlation as 
compared to the standard GW approximation [ 22 ) • The 
cumulant expansion draws from the exact solution of the 
polaron problem [16], and was originally applied to study 
plasmon satellites in core-level spectra m- Importantly, 
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it is also valid in the case of valence electrons, as the 
effects of electron recoil (change of electron momentum) 
upon plasmon emission tend to cancel out [4j. 

In this work we use the formulation of the cumulant 
expansion given by [3] and EL which we will refer to 
as the GW+ Cahk approach. In this formulation only 
the first cumulant is retained in order to describe the 
lineshape of quasi-particles and one-plasmon excitations. 
The case of n-plasmon excitations in > 2) is of little 
interest here since the corresponding spectral signatures 
are damped by the Lang-Firsov factor a n /n\, with a the 
average number of plasmons around the hole jl3j . The 
GW+ Cahk spectral function can be expressed as Eh 

A(k, W ) = £ [A% p (k, U ) + (k, W ) (k, u,)]. (1) 


Here denotes the quasi-particle contribution to the 
Go Wo spectral function: 
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where S' (£") indicates the real (imaginary) part of the 
Go Wo self-energy (2D S3], and e nk the Kohn-Sham eigen¬ 
value. In Eq. ([2]) it is assumed that the off-diagonal ele¬ 
ments of the self-energy are small and can be neglected, 
as is typically the case [23, 331. The term A ( r { in Eq. 0 
is defined as Eh 



H°(k,w) 


/^nk(k-0 /^nk(^nk) (^ 


£-nk) 


d/3 n k 

duj 


(u - £„k) 2 


Enk 


(3) 

where /3„ k (w) = 7r 1 E" k (w)0(p, — to), fi being the chem¬ 
ical potential. This term accounts for interactions be¬ 
tween the photo-hole and one-plasmon excitations [32], 
and its contribution to the spectral function in Eq. 0 is 
to be identified with plasmonic polarons. 

Using Eqs. 0-©. we now investigate the signatures 
of plasmonic polarons in silicon. In Fig. 0a) and (b) we 
report the angle-resolved spectral function of silicon ob¬ 
tained from the Sternheimer-GoWo approach (SGW) and 
SGW plus cumulant (SGW+Cahk), respectively. De¬ 
tails on SGW are provided in the Supplemental Material 
[Mj and Refs. [551155] , For silicon, the experimental (in¬ 
tegrated) photoemission spectrum shown in Fig. 0c) is 
characterized by a broadened plasmonic resonance cover¬ 
ing approximately the energy range from 16 eV to 30 eV 
below the Fermi energy. In this resonance we can clearly 
identify three distinct structures (dashed blue lines). The 
bright energy bands visible in Fig. 0a) for binding ener¬ 
gies between 0-12 eV correspond to the standard quasi¬ 
particle peaks. These peaks result from photoionization 
processes occurring in absence of plasmon excitations, 
and define the ordinary valence band structure of silicon. 
In addition to the quasi-particle features, the spectral 



Figure 1. Angle-resolved spectral function of silicon on a log¬ 
arithmic scale for wavevectors along the T-X high-symmetry 
line, evaluated using (a) the Sternheimer-GW method (SGW) 
and (b) the SGW plus cumulant (SGW+C) approach, (c) 
Measured X-ray photoemission spectrum of silicon (XPS) 
from Ref. JT8] . The blue dashed lines indicate the three fea¬ 
tures discussed in the main text. 


function exhibits a rich structure at binding energies be¬ 
tween 15-30 eV. These structures can be identified with 
plasmonic polarons. These features are present already 
at the G 0 W 0 level, however, their energy range is largely 
overestimated and, thus, not compatible with the plas¬ 
monic features observed in XPS. The inclusion of the cu¬ 
mulant correction in SGW+Cahk moves the plasmonic 
polaron resonances to lower binding energy, improving 
the agreement with the experimental spectrum signifi¬ 
cantly. This improvement reflects the inclusion of higher- 
order exchange-correlation diagrams in the GW+Cahk 
Green function, as discussed in more detail in the Sup¬ 
plemental Material [Ml . 

Unexpectedly, plasmonic polarons exhibit dispersion 
relations which follow closely the ordinary band struc¬ 
ture resulting from the quasi-particle peaks. The strik¬ 
ing similarity between the dispersion of the valence bands 
and the plasmonic structures suggests that we are look¬ 
ing at band structures of plasmonic polarons. Plasmonic 
polaron bands appear as blue-shifted replicas of the or¬ 
dinary valence bands, but they are considerably broader 
and less intense. The comparison of Fig. 0b) and (c) 
suggests that the plasmon satellite of silicon [T8] results 
from the momentum-average of plasmonic polaron bands 
over the first Brillouin zone. For quasiparticles, it is well 
known that the density of states is characterized by singu¬ 
larities (known as Van Hove singularities) at the energies 
for which the first momentum derivative of the quasi¬ 
particle bands vanishes (V k £k = 0). Correspondingly, 
peaks in the density of states can be associated with ex¬ 
tremal points of quasiparticle bands. Figures 0b) and 
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Figure 2. (a) Complete angle-resolved spectral function of monolayer M 0 S 2 on a logarithmic scale for wavevectors along 

the Y-M high-symmetry line, evaluated using the GTT+Cahk approach, (b), (d) Contributions of plasmonic polarons to the 
spectral function in (a): 7r-plasmonic polarons (b) and {ir + er)-plasmonic polarons (d). (c) Quasi-particle band structure of 
monolayer M 0 S 2 extracted from (a) (red solid lines), band structure of 7r-plasmonic polarons extracted from (b) (blue solid 
line), and band structure of (n + <r)-plasmonic polarons from (d) (yellow solid line). Energies are referenced to the valence-band 
top. 


(c) indicate that Van Hove singularities may also have 
a plasmonic origin. In particular, the two experimental 
peaks at 20.4 eV and 24.5 eV and the shoulder at 27.3 eV 
[labelled as 1-3 in Fig.JlJc)] can be attributed to the van¬ 
ishing of the first momentum derivative of the plasmonic 
polaron bands in Fig.Qb). 

Owing to the large plasmon energy, the experimental 
observation of plasmonic polaron bands in silicon may 
be hindered by the low resolution of ARPES measure¬ 
ments well below the Fermi energy. In order to iden¬ 
tify materials in which such polaron bands may be ob¬ 
served, in the following we focus on group IV transition 
metal dichalcogenides (M 0 S 2 , WS 2 , MoSe 2 , and WSe 2 ). 
The EELS spectra of the three-dimensional parent com¬ 
pounds exhibit two distinct features around ~8 eV and 
~22 eV, corresponding to the excitation of it and ir + a 
plasmons, respectively 133 ESI- Since the width of the 
bands arising from transition metal d states and chalco- 
gen p states in these compounds is approximately 7 eV 
[39], possible plasmonic polarons are expected to appear 
between ~8-15 eV. In this energy range the deep S-3s or 
Se-4s bands, located between 12-15 eV [59], dominate the 
spectral function, thereby hindering the identification of 
plasmonic polarons in ARPES also in this case. 

At variance with this scenario, in the case of mono- 
layer TMDs both experimental ||D] and theoretical [41] 
studies reported plasmonic peaks in the EELS spectra 
which are strongly red-shifted with respect to their bulk 
counterpart. For example, in the case of M 0 S 2 and WS 2 
monolayers the 7 r plasmons are found at energies around 
6 eV. We thus expect to observe plasmonic polarons at 
binding energies between 6-13 eV. Since this energy win¬ 


dow matches the band structure gap between the metal- 
d/chalcogen-p bands and the chalcogen s bands, such pla- 
monic polarons should be distinctly observable. 

To examine this possibility on quantitative grounds, 
we calculate the GW +Cahk spectral functions of TMDs 
from first principles [42 j. In order to contain the compu¬ 
tational cost we describe the screening by introducing a 
two-pole approximation for the inverse dielectric matrix, 
as discussed in the Supplemental Material ;5J|. 

Figure [ 2 J a) shows the complete GW+ Cahk spectral 
function of monolayer M 0 S 2 evaluated along the T-M 
high-symmetry line. At binding energies between 0-5 eV 
and 10-15 eV, the spectral function of monolayer M 0 S 2 
exhibits the standard quasi-particle peaks. As in the 
case of silicon, in monolayer M 0 S 2 plasmonic polarons in¬ 
troduce new spectral features in a binding energy range 
where quasi-particle states are absent. To characterize 
these new features, we analyze their energy vs. wavevec- 
tor dispersions using a Lorentzian decomposition of the 
energy profiles. This analysis allows us to disentangle 
the contributions of plasmonic polarons from the quasi¬ 
particle excitations. Figures [2jb) and (d) show the plas¬ 
monic polarons corresponding to the emission of a photo¬ 
electron and the simultaneous excitation of a 7r or 7r + a 
plasmon, respectively. As for silicon, plasmonic polarons 
exhibit a clear energy-momentum dispersion relation that 
leads to the formation of blue-shifted replica of the va¬ 
lence band structure. In particular we find two replicas 
of the valence bands of monolayer M 0 S 2 , one associated 
with the 7T plasmon centered around the binding energy 
~8 eV, and another one associated with the 7r + a plas¬ 
mon around ~19 eV. 
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Owing to the approximate treatment of electron re¬ 
coil effects in the cumulant approach [4, 18], the actual 
broadening might be even larger than in the present cal¬ 
culations. In order to understand whether plasmonic po- 
laron bands could be observed in ARPES experiments it 
is therefore essential to quantify their spectral weight. 
A reliable measure of the intensity of the plasmonic 
bands can be obtained from the average number a„k of 
plasmons emitted during the photoemission process [32| : 
«nk = / w _2 /3„k(w)dw. In the case of the high-lying 
Mo-4d/S-3p bands we obtain a„k in the range 0.08-0.12 
for momenta along the T-M line, whereas a n k is found 
in the range 0.04-0.06 for the lower-lying S-3s bands. 
These results indicate that the formation of plasmonic 
polarons provides an important dissipation channel for 
the ARPES photo-current. In particular, in the case of 
Mo-4d z 2 states at the top of the valence band (the most 
important for electron transport in p-doped M 0 S 2 ) every 
photo-hole is accompanied by 0.13 plasmons per unit cell. 
These estimates are confirmed by a direct integration of 
the spectral function, which shows that plasmonic po¬ 
larons carry on average 12% of the total spectral weight 
of the valence electronic states. For the valence states 
the intensity ratio between the plasmonic polaron band 
and the quasi-particle peak is ~ 0.01-0.02. In addition, 
by considering the metal-d/chalcogen-p bands altogether 
we find that, in the case of monolayer M 0 S 2 , it should be 
possible to extract up to 1.2-10 14 electrons/cm 2 (0.1 elec¬ 
trons/cell) at binding energies inside the band structure 
gap between Mo-4d/S-3p bands and S-3s bands. 

At variance with the case of monolayer M 0 S 2 , for 
the bulk compound our calculations indicate that the 
plasmonic bands overlap substantially with the deep S- 
s states (Fig. [3]). The overlap between ordinary quasi¬ 
particle bands and plasmonic polaron bands should make 
the experimental detection of these new features in bulk 
M 0 S 2 much more challenging than in the case of a mono- 
layer. 

The above results suggest that two key conditions need 
to be satisfied for plasmonic polaron band structures to 
be clearly observed in ARPES: (i) the existence of low- 
energy plasmon excitations in the EELS spectrum, and 
(ii) the presence of a band gap in the valence band mani¬ 
fold. The condition (ii) also contributes to minimize spec¬ 
tral broadening arising from electron-phonon scattering. 

Given these ‘design rules’ it is natural to ask whether 
there exist ‘optimal’ TMDs for observing plasmonic po¬ 
laron bands. To answer this question, we repeated our 
GIE-I-Cahk calculations for the related monolayer com¬ 
pounds MoSe 2 , WS 2 , and WSe 2 - Figure [3] shows that 
plasmonic polaron bands associated with 7 r plasmons 
fit inside the band gap betwen the metal-d/chalcogen- 
p bands and the chalcogen-s bands for all these com¬ 
pounds. In particular we find that monolayer MoSe 2 
should provide an ideal testbench for identifying plas- 
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Figure 3. Energy range of the quasi-particle bands 

(blue boxes) and the 7r-plasmonic polaron bands in mono¬ 
layer M 0 S 2 , MoSe 2 , WS 2 , and WSe 2 , evaluated within the 
GIT+Cahk approach (yellow boxes). Bulk M 0 S 2 is included 
for comparison (MoS 2 [b]). The topmost bands arise from the 
hybridization of transition metal d states and the chalcogenide 
p states (TM-d/Ch-p); the low-energy ones from the chalco¬ 
genide s states (Ch-s). The plasma energies are given in 
eV. 

monic polarons, since in this case the plasmonic struc¬ 
tures exhibit essentially no overlap with the ordinary va¬ 
lence bands (Fig. |3|. 

Since the plasmon energy scales with square root of 
the static dielectric constant e 0 , it should also be pos¬ 
sible to realize plasmonic band structure engineering in 
TMD monolayers by modifying their dielectric screening 
properties. For example this could be achieved by us¬ 
ing different substrates, by building van der Waals het¬ 
erostructures of TMDs [33] , using doping [33] , or mechan¬ 
ical deformation [33] [33]. All these modifications would 
leave the ordinary valence band structure essentially un¬ 
affected (on the eV scale), while producing significant 
shifts in the binding energy of plasmonic polarons. 

In conclusion, using state-of-the-art first-principles 
GW plus cumulant calculations, we have shown that 
electron-plasmon coupling leads to the formation of plas¬ 
monic polaron band structures, and that group IV tran¬ 
sition metal diclialcogenide monolayers (in particular 
monolayer MoSe 2 ) should provide a unique opportu¬ 
nity for observing these new features in ARPES exper¬ 
iments. Remarkably, we have found that these plas¬ 
monic band structures exhibit dispersion relations which 
closely follow the ordinary valence bands and, similarly 
to quasiparticle bands, lead to the formation of Van Hove 
singularities in the density of states. More generally, 
plasmonic-polaron bands emerge as novel spectroscopic 
signatures of electron-plasmon coupling, which may con¬ 
tribute to unravel the complexity of ARPES measure¬ 
ments. The emergence of plasmonic polarons is not lim¬ 
ited to semiconductors, and may also prove useful to 
rationalize the electronic structure of materials charac¬ 
terized by well defined plasmonic excitations, including 
for example metals, as well as d- and /-electron systems 
mm- In fact, our work raises the question on whether 
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the concept of plasmonic polarons may be systematically 
employed in the interpretation of ARPES spectra of com¬ 
plex systems. As an example, while our work focused on 
the simplest case of ‘isolated’ plasmonic polaron bands, 
the crossing of plasmonic bands and high-lying valence 
bands may reveal the emergence of band-branching ef¬ 
fects, in analogy with the polaron problem in electron- 
phonon physics [49]. 

This work was supported by the Leverhulme Trust 
(Grant RL-2012-001) and the European Research Coun¬ 
cil (EU FP7 / ERC grant no. 239578 and EU FP7/grant 
no. 604391 Graphene Flagship). Calculations were per¬ 
formed at the Oxford Supercomputing Centre and at the 
Oxford Materials Modelling Laboratory. 
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